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Abstract. We provide evidence in support of a recent proposal by Astrakharchik 
et al for the existence of a super Tonks- Girardeau gas-like state in the attractive 
interaction regime of quasi-one-dimensional Bose gases. We show that the super Tonks- 
Giradeau gas-like state corresponds to a highly-excited Bethe state in the integrable 
interacting Bose gas for which the bosons acquire hard-core behaviour. The gas-like 
state properties vary smoothly throughout a wide range from strong repulsion to strong 
attraction. There is an additional stable gas-like phase in this regime in which the 
bosons form two-body bound states behaving like hard-core bosons. 
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A sequence of striking advances has seen the experimental realization of quantum 
gases of ultracold atoms in optical lattices. Experiments in trapping one-dimensional 
(ID) Bose gases have been remarkably successful [TJ El El HI Ej • The effective ID coupling 
constant #id can be tuned essentially from the weak coupling regime to the strong 
coupling regime via anisotropic confinement or Feshbach resonance. The ratio of the 
average interaction energy to the kinetic energy, 7 = m A 1D , is used to characterize the 
different physical regimes of the ID quantum gas. Here m is the atomic mass and n is the 
boson number density. For the regime in which the interactions are weakly repulsive, i.e., 
7< 1, the ID degenerate Bose gas acts like a quasi Bose-Einstein condensate. As the 
strength of the repulsive interaction tends to infinity the bosons behave like impenetrable 
fermions and the system is known as the Tonks- Girardeau gas. In this regime, the ground 
state properties of the ID homogeneous Bose gas can be investigated using the wave 
function of noninteracting fermions [H] as a special limit of the Bethe Ansatz solvable 
Bose gas [7j. On the other hand, observation of the formation of bright solitons in quasi- 
1D trapped condensates [HI E] has renewed interest in the exactly solved ID model of 
interacting bosons in the attractive regime fUl ■ It was shown that attractive bosons 
confined to a finite system with periodic boundary conditions can undergo a transition 
from uniform condensation to bright solitons (TTJ [T2 [TH] . These studies reveal subtleties 
in the signature of bosons in ID many-body systems [El [121 IH3 U3 UH EH EH! 

It is very interesting to consider what happens to the ID Bose gas when the coupling 
strength is abruptly changed from strong repulsive interaction to strong attractive 
interaction - does the gas suddenly collapse or does it inherit a stable gas-like phase 
from the hard-core bosons? A recent proposal |2*T) 122] by the group at BEC-INFM in 
Italy suggests the existence of a super Tonks- Girardeau gas for the strongly attractive 
regime. They argued that this new gas-like state can be stable in a wide range of 
strongly attractive interaction strength. In using a variational Monte Carlo calculations 
they found that this gas-like state can be described by a gas of hard rods with strong 
attraction, which yields a positive compressibility. In this regime the ID Bose gas is 
more strongly correlated than the Tonks-Girardau gas - they thus referred to it as the 
super Tonks-Girardau gas. Here we show that this super Tonks-Girardau gas-like state 
corresponds to a highly-excited state in the Bethe Ansatz solution of the interacting Bose 
gas for strongly attractive interaction. We also predict the existence of another stable 
gas-like phase where the bosons form two-body bound states with large binding energy. 
The existence of these gas-like states is due to the fact that a large Fermi-pressure-like 
kinetic energy inherited from the Tonks- Girardeau gas cannot be quenched instantly 
when the interaction is switched from strong repulsion to strong attraction. This is 
analogous to a system in classical mechanics where a moving ball is able to reside on 
a step of a declining staircase in a potential well if the friction between the ball and 
the surface of the steps is strong enough. In addition, for weakly attractive interaction 
the bosons form a many-body bound state which leads to the ground state energy per 
particle being proportional to the linear density n. 
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The interacting <5-function Bose gas [2H] is described by the Hamiltonian 

h 2 N d 2 

i=l 1 l<i<j<N 

in which iV particles are constrained by periodic boundary conditions on a line of length 
L. Here #id = h 2 c/m is an effective ID coupling constant with scattering strength 
c. If the gas is harmonically trapped along the axial direction, the scattering strength 
is determined by c = 2/|oid|, where aiD is the effective ID scattering length. The 
eigenvalues are given by E = ^ Yuf=i kj > where the quasi-momenta kj satisfy the 
Bethe equations 

exp(i%L) = - n g ^ + 1C , for j = l,...,iV. (2) 
^-^ kj — kp — ic 

The Bethe roots fcj are known to be real for repulsive interactions, c > 0, where the 
model has been extensively studied (see also, Refs. [2H 123 HH 123 12H1 121] and references 
therein). A significant hallmark of ID bosons is that the quasi- momentum distribution 
satisfies a universal semi-circle law for relatively weak interaction strength 0EO1I2Z]- As 
the interaction strength increases the distribution function becomes flat and the bosons 
behave like impenetrable free fermions. More generally, the long-wavelength properties 
of the ID Bose gas can be described by a Luttinger liquid due to the universality in the 
low-energy excitations, i.e., gapless excitations with a linear spectrum in the low-energy 
excitations and power-law decay in the correlations [SU E21 EH] ■ 

On the other hand, the situation is far from completely understood for attractive 
interactions, c < 0, where the roots form bound states in the complex plain. It is 
believed that for attractive interaction there is no saturation due to the absence of a 
hard core in the strong coupling limit. McGuire [21] argued that the ground state energy 



for iV particles in infinite space is Eq = —j^c 2 N(N 2 — 1), which agrees with the result 
later obtained from the string solution [23] . This implies that there is no thermodynamic 
limit. However, Astrakharchik et al |21[l2*2*] have suggested the existence of a new stable 
gas-like phase in the ID Bose gas for strongly attractive interaction. Using variational 
Monte Carlo calculations they found that for a certain range of strong coupling, the 
energy per particle coincides with the energy 

N 6m V c J K ' 

of a gas of hard rods. This leads to a positive compressibility. In this regime there exists 
an energy barrier which separates this gas-like state from cluster-like bound states [33J. 
It was also shown that the short range correlations for the super Tonks-Girardeau gas- 
like phase are much stronger than in the usual Tonks-Girardeau gas. 

Before discussing the strongly attractive regime, consider first the weakly attractive 
case. Here we see that the exactly solved ID model of attractive interacting bosons 
does saturate and the ground state energy for the gas phase changes smoothly as the 
interaction varies from weakly repulsive to weakly attractive. Following [2Z] , we find that 
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the roots of the Bethe equations for weakly attractive bosons satisfy a set of algebraic 
BCS-like equations of the form 



where the summation excludes j = £ and dj = 0, ±1, ±2, . . . denote excited states. 
The ground state has zero total momentum, with dj = for j = 1, . . . , N, with the 
roots forming an iV-body bound state. These equations are the same as for the weakly 
repulsive regime (HOI 1211, although now the roots of equation (JU) are all purely imaginary 
for the ground state as c < 0. They are given by Hermite polynomials of degree N, 
i.e., iJAr(y^fc) = |S0I - Alternatively, the roots may be given in terms of Laguerre 
polynomials [21]. The ground state energy per particle is E /N = c(N — 1)/L, which 
follows directly from (J1J) (see also Ref. [S3)- I n this way the ground state energy per 
particle in the weakly attractive and repulsive regimes is seen to pass continuously 
through the origin at c = 0. Since our analysis is based on finite systems, it avoids 
the subtlety of a singular integration kernel in the Lieb-Liniger integral equation in the 
limit c — * 0. Our findings agree with the previously reported numerical results for the 
attractive regime [3H ESI EEO • 

An alternative approximation of the Bethe equations for weak interaction consists in 
observing that in this regime, each configuration of the Bethe-numbers can be described 
by N c clouds of quasi-momenta, centred around k^ = 2-nj/L, where integer j / 
corresponds to excited states. The momenta of each cloud lie parallel to the imaginary 
(on the real) axis for c < (c > 0). More quantitatively, let us label the fc's by two 
indices j and u, where j = l,...,N c and v counts the particles rij in each cloud j, 
v = 1, . . . , rij. Then 



as the "centre of mass" of the j-th cloud in quasi-momentum space and (3l describes 
the position of the individual quasimomenta within one cloud around its centre of mass. 
By linearizing the Bethe equations including 0(c) terms, one obtains 



Note that according to (|7|) . the momenta within one cloud are given by roots of Hermite 
polynomials, centred around the cloud centre of mass. Only this centre of mass is 
influenced by the presence of other clouds, cf. ©. Thus the linearization of the Bethe 
equations is analogous to a dipole expansion in electrostatics, with charges placed at the 
position of the quasimomenta. Note that Eqs. (JBJ) and (JZj) coincide with Eq. (J1J) only for 
the ground state. The reason is that for excited states, Eq. (JU) contains terms of order 







(6) 



(7) 
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higher than 0(c). It turns out that Eq. (jlj) yields a better approximation for c < 0, 
whereas Eqs. (JHJ) and (JJJ) are better suited for c > 0. It is straightforward to calculate 
the energy E c i ou d from Eqs. (jBJ) and (JZJ), 

E doud = [pf] n i + x N ( 2N - 1) - E 4 ( 8 ) 

As an example, consider iV = 12 weakly interacting bosons. Tables 1 and 2 show 
a direct comparison between the ground state energy and the leading eight excitations 
obtained from the various equations. These are the energies obtained directly from 
numerical solution of the Bethe equations (0), the BCS-like equations and the cloud 
equations (JHJ)-©- Also shown are the energies obtained from the direct diagonalization 
[TT] of the truncated Hamiltonian. The overall agreement is excellent, with even very 
close energy values clearly distinguishable by their characteristic root distributions and 
total momentum. 

The momentum density distribution for the ground state is along the imaginary 
axis, i.e., the Bethe roots are of the form ikj for j = 1, . . . , N, and given by the semi- 



circle law n(k) = ^=yl — 4 ^ n2 ■ This is analogous to the behaviour in the weakly 
repulsive case [7J|3UJ|2I]. The momentum distribution also becomes flat if the interaction 
strength is increased. The stronger the interaction strength, the larger the momentum 
distribution region. This implies a signature for the ID Bose gas in the weak attractive 
regime 7 <C 1/N. 

We turn now to strongly attractive interaction. From the exact analytic result for 
two bosons with attractive delta interaction, if c < — 2, the roots for the ground state 
tend to i|. Without loss of generality, we focus on N = 2M bosons with Lc <C — 1 for 
even M. Then for general N, the roots of the Bethe equations (J2J) for the ground state are 
approximately given 

by |23 123 k±j « ±i (f (2M - 2j + 1) + Sj) , j = 1, . . . , M, where 
Sj is small and negligible for Lc <C — 1. In this case the wave function is approximately 
given by [T^] 



^(xi,...,Xiv) ~ A/"exp< - \xj-Xi\y. (9) 




yj (n,— 1)1 (n-1) 

Here N = — - c is a normalization constant. The eigenvalue of the Hamiltonian 



for this state is Eq = —j^c 2 N(N 2 — 1), which is McGuire's soliton-like ground state 
energy [23] . 

On the other hand, if the interaction strength is abruptly changed from strongly 
repulsive to strongly attractive, the gas-like phase may vary smoothly against the cluster- 
like phase |2U EH]- We see this in the context of the exactly solved model from a 
highly-excited state in which all Bethe roots are real and symmetric about the origin, 

with ±A; 2m _i, m = 1, . . . , N/2, where hi ~ ^ ^1 + . These roots can be directly 
obtained from the asymptotic expansion of the Bethe equations (J2J) for L\c\ 3> 1. The 
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energy of this state follows directly from this asymptotic solution, with 

E ft 2 I, „ ,n 2 f 2n\~ 2 . . 

^^^(N 2 -l)- 1 + - (10) 



N 2m 3 v ' L 2 V c / 

which coincides with the result ((5J) for the gas of hard- rods j^U EH] • Here the condition 
■jrj^ <C 1 is required. For harmonic trapping along the axial and the radial directions, 
this condition becomes N <C L/am, where L is proportional to the characteristic 
length a z = ^h~/muj Z) i.e. L oc \fNa z . Therefore the strong coupling condition is 
equivalent to a restriction, iV <C a 2 /a 2 D , on the number of bosons, which is similar to 
the requirement for the formation of bright solitons in attractive condensates [HIE]- Here 
oj z is the trapping frequency along the axial direction. We obtain the sound velocity 
v s =vf/(1 + -) 2 which is positive for strongly attractive coupling. This indicates that 
the compressibility [2T] of the system is positive, corresponding to a stable gas-like 
phase. The Fermi velocity is Vp = The Luttinger liquid parameter characterizing 
the density correlation rate is K = (1 + -) 2 , indicating an enhancement of the density- 
density correlations in this super Tonks- Girardeau gas. For example, the one-particle 
correlations g(x) = ( , ip^(x)ip(0)) oc I/x^k. In order to see the enhancement of the 
correlation function in the super Tonks-Girardeau phase, we write K = (1 — |aiD |n) 2 , 
which is smaller than in the Tonks-Girardeau phase, where K — (1 + |aio |^) 2 - 

The local two-body correlation function = Aix 2 n 2 /37 2 (1 + -) 3 measures the 
probability of collisions between two particles. The interaction energy of the ultracold 
bosons depends on coherent collisions. In the strong coupling limit, the probability of 
collisions is very small due to decoherence of the wave functions, i.e., g^ — > 0. Therefore 
the interaction weakly affects the energy per particle in the strong coupling regime due 
to the existence of kinetic pressure. Figure 1 shows a plot of the energy per particle (fTU|) 
and the velocity versus the gas parameter na 1D . We see that they vary smoothly as the 
interaction varies from strongly repulsive to strongly attractive. The existence of Fermi- 
pressure-like kinetic energy is the main reason for a stable gas-like state throughout the 
whole of the crossover from strong repulsion to strong attraction. 

Nevertheless, there is another stable gas-like state where the Bethe roots form 
two-body bound states, namely with complex conjugate pairs of the form ±&2 m -i ~ 
a m ± i (| + S m ) , m = 1, . . . , with 

(2m- 1)tt / 5MV 1 , , 



(12) 



2L \ 2L C/ 

_ 2 M-2 c M e §L 



< nS 2 IM) 2 - «) 2 I + 2^- 2 Lc^' 
The case I — m is excluded in the product in p2jl . The energy per particle follows as 

E J ^ 2 _J_c 2 _ tt 2 (M 2 - 1) 2 M c M+1 e§ L (l + f^)^ 1 1 
N*2m~{ 4 + 12L 2 (l + |M) 2 ^rT(M) /' (13) 

Here T denotes the gamma function. In the RHS of the above equation the first term 
is the pair binding energy. The second term contains the kinetic energy and interaction 
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energy. The last term is the correction from S m in (J 12)) . This correction term is negligible 
in the thermodynamic limit. Thus we have 

-2' 




(14) 

This state is an analogue of the fully-bound fermion pairs in the ID Fermi gas where 
the bound pairs cannot condense due to Fermi pressure [23 EH1 HUH HOI SJj • Those pairs 
can be thought of as hard-core bosons. For this state the sound velocity is effectively 
given by v s = \vp(l + ^)~ 2 > 0, implying the stability of this gas-like state. The 
Luttinger liquid parameter is K — 4(1 + ^) 2 . The enhancement of K in this phase 
leads to a very slow decay of the phase correlations jSDE21, but less strongly correlated 
than the super Tonks-Girardeau gas. The sound velocity decreases with the formation 
of multi-body bound states. The vanishing of the sound velocity implies that the system 
is unstable in multi-soliton-like states where the system lacks hard-core behaviour. But 
some soliton-like states which lead to positive compressibility would be stable. This is 
a manifestation of the formation of matter-wave soliton trains in quasi-lD attractive 
condensates [HJ Ej • 

We conclude that some gas-like states could be stable against cluster-like states 
if the interaction strength abruptly changes from strong repulsion to strong attraction. 
The gas can inherit hard-core behaviour due to the existence of large Fermi-pressure-like 
kinetic energy from the strongly repulsive interaction. This regime could potentially be 
reached in experiments via anisotropic confinement in which a resonance of the effective 
ID scattering length takes place. 
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Table 1. The ground state energy per particle (j = 0) and the first nine excitations 
(j = 1, . . . , 9) for N = 12 weakly interacting bosons with c = —0.05 and L = 6. Also 
shown is the total momentum P of each state. E e ^ ct denotes the energies obtained 
directly from numerical solution of the Bethe equations (j2J). E BCS denotes energies 
obtained from solutions of the BCS-like equations (@J. E cloud denotes energies obtained 
from solutions of the cloud equations I©-©- For further comparison, E nurn , denotes 
energies obtained from direct diagonalization of the truncated Hamiltonian. 



j 


P 


rp 

'exact 


E B cs 


-'-''cloud 


rp 

■'-'num. 








-0.11028 


-0.11000 


-0.10000 


-0.11021 


1 


tt/3 


0.96786 


0.96809 


0.96829 


0.96794 


2 





2.04772 


2.04791 


2.04825 


2.04782 


3 


2tt/3 


2.04936 


2.04955 


2.04991 


2.04946 


4 


tt/3 


3.13094 


3.13109 


3.13153 


3.13433 


5 


71 


3.13422 


3.13438 


3.13487 


3.13432 


6 





4.21587 


4.21599 


4.21649 


4.21599 


7 


2tt/3 


4.21751 


4.21764 


4.21816 


4.21763 


8 


4tt/3 


4.22244 


4.22260 


4.22316 


4.22256 


9 


2vr/3 


4.25788 


4.25810 


4.25816 


4.25797 



Table 2. Same as for Table 1, but with L = 24. 



j 


P 


rp 

-^cxact 


E B cs 


rp 

cloud 


rp 

■'-'num. 








-0.02779 


-0.02750 


-0.02750 


-0.02772 


1 


tt/12 


0.03601 


0.03624 


0.03646 


0.03609 


2 





0.10029 


0.10047 


0.10083 


0.10039 


3 


7r/6 


0.10067 


0.10086 


0.10125 


0.10077 


4 


tt/12 


0.16542 


0.16557 


0.16603 


0.16554 


5 


tt/4 


0.16619 


0.16635 


0.16687 


0.16631 


6 





0.23102 


0.23114 


0.23166 


0.23115 


7 


7r/6 


0.23141 


0.23153 


0.23207 


0.23154 


8 


tt/3 


0.23258 


0.23271 


0.23332 


0.23270 


9 


7r/6 


0.24178 


0.24202 


0.24207 


0.24188 
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Figure 1. The energy per particle in units of H 2 n 2 /2m versus the gas parameter 
naiD for strong attractive and strong repulsive interaction. Here gin = 2H 2 /man). 
The inset shows the sound velocity versus the gas parameter. The gas-like state varies 
smoothly due to the existence of Fermi-pressure-like kinetic energy if the interaction 
strength is abruptly changed from strongly repulsive to strongly attractive. 



